Size segregation in granular media induced by phase transition 
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In order to study analytically the nature of the size segregation in granular mixtures, we introduce 
a mean field theory in the framework of a statistical mechanics approach, based on Edwards' original 
ideas. For simplicity we apply the theory to a lattice model for hard sphere binary mixture under 
gravity, and we find a new purely thermodynamic mechanism which gives rise to the size segregation 
phenomenon. By varying the number of small grains and the mass ratio, we find a crossover from 
Brasil nut to reverse Brasil nut effect, which becomes a true phase transition when the number of 
small grains is larger then a critical value. We suggest that this transition is induced by the effective 
attraction between large grains due to the presence of small ones (depletion force). Finally the 
theoretical results are confirmed by numerical simulations of the 3d system under taps. 

PACS numbers: 45.70.Mg, 64.75.+g, 05.50.+q 
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In the last decades a great attention has been devoted 
to the study of the problem of vertically shaken granu- 
lar mixtures under gravity. It was observed that such 
systems depending on the control parameters can mix or 
segregate their components spontaneously according to 
a mechanism which is still largely unclear, although of 
deep practical and conceptual relevance P, 0. Rosato 
et al. showed that large grains surrounded by a sea 
of smaller ones rise to the top when subjected to verti- 
cal shaking. This is the well known "Brazil Nut Effect" 
(BNE), while the opposite one (i.e. large grains on the 
bottom and small ones on the top) is known as "Reverse 
Brazil Nut Effect" (RBNE) HQ. The BNE was origi- 
nally explained [3j in terms of a geometric effect, based 
on "percolation" arguments. Along with geometry, dy- 
namical effects, such as convection jg] or inertia @, were 
shown to play a crucial role. The hydrodynamic equa- 
tions for binary mixtures Q were also applied, and in 
this framework the regimes of BNE and RBNE Q were 
found to depend on the presence of inelastic dissipation. 

Recent results finally outlined that segregation pro- 
cesses may involve global mechanisms, such as condensa- 
tion p| or phase separation 01 ■ This suggested a change 
of perspective on the issue and the idea to formulate 
a statistical mechanics description of these phenomena 

II IH 113 To this aim > in Ref - S II EH the interplay 
between size and mass was studied and a phase diagram 
for BNE/RBNE transition based on the competition be- 
tween percolation and condensation was proposed. 

In the present paper, applying a statistical mechan- 
ics approach for non-thermal systems |l'il If4| . a new 
purely thermodynamic mechanism is found responsible 
for size segregation. We follow a statistical mechanics 
approach for granular materials under taps, which was 
recently developed [lij on the basis of Edwards' original 
ideas 0, Q|. This approach postulates that time av- 
erages coincide with suitable ensemble averages over the 
mechanically stable states (i.e. those where the system 



is found at rest). Such hypothesis was shown 0] to hold 
with good approximation in a lattice model for granular 
mixtures under taps. In order to study size segregation, 
in this framework we apply a mean field theory to such 
model. 

We find that, by fixing the number of grains of the 
two species, it is possible to obtain a smooth crossover 
from BNE to RBNE by changing the mass ratio of the 
grains. However if the number of small grains increases 
above a critical value, a sharp transition occurs resulting 
in a phase separation between the two species: Verti- 
cal size segregation is induced by the presence of gravity 
which forces the heavier phase to move downward. At 
the critical point, critical fluctuations are found typical 
of a second order phase transition. We suggest that this 
demixing phase separation is due to the presence of the 
depletion force. This is an effective attraction between 
large grains, obtained by tracing over all the degrees of 
freedom of small grains. The strength of the interaction 
is proportional to the number of small grains. The deple- 
tion force, typically present in thermal binary mixtures 
[lfij . in fact may be also present in 3d granular mixtures 
under gravity [17] . The theoretical results are confirmed 
by numerical simulations of the 3d system. 

The model - The model E3 is a hard sphere bi- 
nary mixture under gravity made up of two species, 1 
(small) and 2 (large), respectively with diameter ao = 1 
and \/2ao, and masses mi and m^. In order to sim- 
plify the calculations we divide the space in cubic cells 
of linear size ao whose vertices can be or not occupied 
by grains. Each site of the lattice is labeled by the two 
indices i — 1, . . . , L 2 and z = 1, . . . , L z , where L and L z 
are respectively the horizontal and vertical dimensions 
of the lattice. A microstate of the model is completely 
identified by the occupancy variables {hJ}, associated to 
each site of the lattice: n\ = 0, 1, 2 respectively if the site 
(i, z) is empty, filled by a small grain or by a large one. 
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In terms of these variables the model Hamiltonian is 

H = Hhc + migHi + m 2 gH 2 , (1) 

where Hi = 52 iiZ z5 n *i, H 2 = J2i, z z ^np are respec- 
tively the heights of the two species, and Hhc is the hard 
core potential, preventing two nearest neighbor sites to 
be both occupied if at least one contains a large grain. 

In Ref.s hj it was shown that the model, Eq. QJ, 
under taps is described with good approximation by a 
statistical mechanics approach based on Edwards' origi- 
nally ideas. In particular it was found that the weight of 
a microscopic state, {nf}, is given by 

p/S n z \^j = n({n z j) • e - H Hc({ni})~l3im 1 gH 1 -p 2 m2gH 2 

(2) 

where II({n|}) is a projector over the "mechanically sta- 
ble" states, namely n({nf}) = 1 if {nf} is mechanically 
stable, and II({nf}) = otherwise, f3\ and j3 2 are the 
variables conjugated respectively to the gravitational en- 
ergies of the two species, and the two "configurational 
temperatures" T\ — /3-j -1 and T 2 — (i^ 1 are increasing 
functions of the tap amplitude. The system partition 
function is given by Z = J2{n z \ P{{ n i }) where the sum 
is over all the microstates but, due to the projector, only 
the mechanically stable ones are taken into account. In 
the following calculations we adopt a simple definition of 
"mechanical stability" : a grain is considered stable if it 
has a grain underneath |18j. 

Mean field approximation - In the present paper the 
partition function, Z, is evaluated in mean field approx- 
imation: we consider a generalization of Bethe-Peierls 
method for anisotropic systems (due to gravity) already 
used in previous papers 0, 0, |2jj. In particular, we 
consider a Bethe lattice made up by L z horizontal lay- 
ers (i.e., z € {1, ...,L Z }) whose edges can be occupied by 
grains. Each layer is a random graph of given connectiv- 
ity, k — 1 = 4. Each site in layer z is also connected to its 
homologous site in z — 1 and z + 1 (the total connectivity 
is thus fc + 1 = 6). We use the grancanonical ensemble by 
introducing two chemical potentials, fi\ and fi 2 , relative 
to species 1 and 2. 

In the following we give the main ideas of the cal- 
culations and refer to a longer paper [2(| for the de- 
tails (see also 21] ). The Bethe-Peierls recursion equa- 
tions, which allow to calculate the partition function, 
are written in terms of the local "cavity fields" defined 
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(with n = 1,2). Here Z t 



and Zn's^ 



are the partition 



functions of the "side" branch restricted respectively to 
configurations where the site (i, z) is empty or filled by a 

particle of specie n; analogously, Zn^\ Zq'*^ and Z^'u 

(resp. Z^^\ Zq'^ and Z^'d ) are the partition functions 
of the "up" (resp. "down") branch restricted to configu- 
rations where the site (i, z) is filled by a grain of species 



n, empty with the upper (resp. lower) site empty and 
empty with the upper (resp. lower) site filled by a grain 
of species n. 

The calculations are limited to the "fluid" phase corre- 
sponding to a solution of Bethe-Peierls equations where 
local fields in each layer are site independent j^. Such 
a solution, characterized by horizontal translational in- 
variance, is given by the fixed points of the following 
equations: 
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3 /3„( Mn -m„ S 2)_ The func _ 



tions appearing in Eq. have a simple expression in 
terms of the local fields l23l. From the local fields we 
calculate the free energy [24j, F, and all the quanti- 
ties of interest, such as density profiles of small and 
large grains, cri(z) and cr 2 (z), the grain number per 
unit surface, N n — ^2 z cr n (z), and their average height, 

h n = (Z n ) = J2z Z(T n{z)/J2z a n(z) (with U = 1,2). 
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FIG. 1: (color online). Main frame as a function of S in 
the mean field approximation for fixed value of N2 = 0.685, 
and Ni — 0.5 (blue curve), Ni = 0.9 (green curve), 
Ni = TVic = 1.37 (red curve) and Ni = 2.4 (black curve) 
(TVi increases along the arrow). Inset: The susceptibility, Xi 
as a function of 5 for fixed value of A^2 = 0.685, and Ni = 0.5 
(blue curve), Ni = 0.9 (red curve) and Ni — Ni c = 1-37 
(black curve) (Ni increases from bottom to top). 

In the following instead of using the chemical potential 
variables, fii and /i2, we use the conjugate variables, N% 
and N 2 , respectively the number of the small and large 
grains per unit surface. For simplicity we fix the config- 
urational temperatures, T± = T 2 = 1, and the number 
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of large grains, N%, and vary N%. In order to study the 
BNE/RBNE transition we also vary the masses of the 
two species, mj and 7712 (keeping constant mj + m 2 ), 

This 



— h\ — hi 



hi+h 2 



2m± — ni2 



and measure the height difference, 
quantity, that is thus a function of N-i and 8 — , 
is positive when small grains are in average above large 
grains (RBNE) and negative otherwise (BNE). 
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FIG. 2: (color online). The coexistence curve (black curve on 
the left) and the maximum of the susceptibility (blue curve 
on the right) as a function of N^ 1 . 

In the Main Frame of Fig. [TJ ^(S,Ni) is plotted for 
fixed Ni as a function of 8: The different curves corre- 
spond to increasing values of N\ (along the arrow). We 



Ah 



is a contm- 



find that for low enough values of Ni 
uous monotonous decreasing function of S, due to en- 
tropic and buoyancy effect. By increasing Ni above 
some critical value (Ni c ~ 1.37) Ah. instead exhibits 
first order transition with a finite jump from a positive 
value to a negative one at 5 = 8 CO ex{Ni), where a BNE 
phase and a RBNE one coexist with the same free en- 
ergy. As shown in Fig. El the jump becomes smaller and 
smaller approaching Ni c from above, and goes to zero 
continuously at Ni c (the critical point is (<5 C , Ni c ), where 
5 C = lim Nl ^ Nla Sco^Ni)). 

Along the coexistence curve near the critical point in 
both phases the order parameter 4^ ~ (N^. 1 — N^ 1 ) 13 , 
where, solving the mean field equations numerically, p is 
found consistent with the mean field value (3 — 0.5. 

As in any second order phase transition, the fluc- 
tuations of the order parameter diverge at the criti- 
cal point. In the Inset of Fig. the susceptibility, 
X(8,Ni) = N tot {{Ah 2 ) - (Ah) 2 ) (where N tot is the total 
number of grains) is plotted as a function of 6 (the curves, 
from bottom to top, correspond to increasing values of 
Ni). As Ni approaches Ni c , the susceptibility diverges 
at 8 = 8 C . For Ni < Ni c , x($,Ni) displays a maximum 
(for the value of 8 where Ah is zero), which diverges as a 
power at iVi c , {N± —Ni c )~ J with an exponent 7 found, 



solving numerically the mean field equations, consistent 
with the mean field value 7 = 1.0 (see Fig.[2J. 

In the present paper we have chosen as free parameters 
Ni and 6: 8 plays the role of an external field, and Ni 
plays the role of an inverse temperature, since it controls 
the amplitude of the depletion force 0] . To confirm this 
idea, we have studied the case of mixture of grains with 
the same radii (R2/R1 = 1), where the depletion force is 
completely absent. We find that the critical point disap- 
pears, and a smooth crossover from positive to negative 
values of Ah, varying the control parameters Ni and <5. 

Monte Carlo simulations - The theoretical results ob- 
tained in mean field approximation arc confirmed by 
Monte Carlo simulations of the model, Eq. QJ, under 
taps 0, Hj| on a cubic lattice. In Fig. |3| the height 
difference, and the susceptibility, are plotted as 
function of i5 for fixed N2 and increasing values of Ni. 




FIG. 3: (color online). Main frame 4^ as a function of 5 
in the Monte Carlo simulations on the cubic lattice for fixed 
value of ./V2 = 1.2, and Ni = 0.6 (blue curve and filled cir- 
cles), Ni = 0.8 (green curve and empty circles), Ni = 1.0 
(red curve and stars) and Ni — 1.2 (black curve and trian- 
gles). Inset The susceptibility, x> as a function of <5 for fixed 
value of JV2 = 1.2, and Ni = 0.6 (blue curve and filled cir- 
cles), N\ = 1.0 (red curve and stars) and Ni — 1.2 (black 
curve and triangles) 

As found in the mean field approximation, at fixed 



value of Ni, 



Ah 



as a function of 6 displays a steeper 



and steeper slope (and the susceptibility, Xi a larger and 
larger maximum) as Ni is increased. A power law di- 
vergence of the maximum of the susceptibility is found 
at Nic — 2.2, with an exponent 7 ~ 1.6. As usually 
in mean field approximation, the critical point (N^, 1 ) is 
over-estimated. 

As Ni is further increased, Ah. obtained by varying 8 at 
a given rate, becomes irreversible, and a strong hysteresis 
appears, signaling the presence of a first order transition. 

Conclusions- We focus on the problem of the verti- 
cal size segregation in binary granular mixtures. In the 
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framework of a statistical mechanics approach to granu- 
lar media, based on Edwards' original ideas, we apply a 
mean field theory to a simple hard sphere lattice model 
under gravity, and explain BNE and RBNE with purely 
thermodynamic argument. 

We find that by varying the control parameters a 
transition from BNE to RBNE can occur as a smooth 
crossover or as a sharp transition, depending whether 
one is above or below the critical point. The presence 
of such a critical point is manifested by the divergence 
of the height fluctuations of the two species. Therefore 
we suggest that it might be important to measure in ex- 
periments these fluctuations in the vicinity of the transi- 
tion from BNE to RBNE. The larger the fluctuations the 
closer to the critical point are the control parameters. 

We would like to thank Hans Herrmann for many inter- 
esting discussions and suggestions. Work supported by 
EU Network Number MRTN-CT-2003-504712, MIUR- 
PRIN 2004, MIUR-FIRB 2001, CRdC-AMRA, INFM- 
PCI. 



[1] J. Bridgwater, Chem. Eng. Sci. 50, 4081 (1995). J.M. 
Ottino and D.V. Khakhar, Ann. Rev. Fluid Mech. 32, 
55 (2000). T. Shinbrot and F.J. Muzzio, Physics Today 
53, 25 (2000). Rosato ct al., Chem. Eng. Sci. 57, 265 
(2002). 

] M.E. Mobius, X. Cheng, G.S. Karczmar, S.R. Nagel, 

H.M. Jaeger, cond-mai/0405328. 
.] A. Rosato, K.J. Strandburg, F. Prinz, R.H. Swendsen, 

Phys. Rev. Lett. 58, 1038 (1987). 
] D.C. Hong, PV. Quinn, S. Luding, Phys. Rev. Lett. 86, 
3423 (2001). 

] A. P. J. Breu, H.-M. Ensner, C. A. Kruelle, and I. Re- 

hberg, Phys. Rev. Lett. 90, 014302 (2003) 
.] J.B. Knight, H.M. Jaeger, S.R. Nagel, Phys. Rev. Lett. 

70, 3728 (1993). 
'] T. Shinbrot and F.J. Muzzio, Phys. Rev. Lett. 81, 4365 
(1998). 

,] J.T. Jenkins and D.K. Yoon, Phys. Rev. Lett. 88, 194301 
(2002). 

'] L. Trujillo, M. Alam and H. Herrmann, Europhys. Lett. 

64, 190 (2003). 
i] P.M. Reis and T. Mullin, Phys. Rev. Lett. 89, 244301 
(2002). P. M. Reis, G. Ehrhardt, A. Stephenson and T. 
Mullin Europhys. Lett. 66, 357 (2004). 
] J.A. Both and D.C. Hong, Phys. Rev. Lett. 88, 124301 
(2002). 

] M. Tarzia, A. Fierro, M. Nicodemi, and A. Coniglio, 

Phys. Rev. Lett. 93, 198002 (2004). 
.] S.F. Edwards and R.B.S. Oakeshott, Physica A 157, 
1080 (1989). A. Mehta and S.F. Edwards, Physica A 157, 



[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 
[10 

[11 

[12 

[13 



1091 (1989). 

[14] M. Nicodemi, Phys. Rev. Lett. 82, 3734 (1999). A. Bar- 
rat, J. Kurchan, V. Loreto and M. Sellitto, Phys. Rev. 
Lett. 85, 5034 (2000). J.J. Brey, A. Prados, B. Sanchez- 
Rey, Physica A 275, 310 (2000). A. Coniglio and M. 
Nicodemi, Physica A 296, 451 (2001). D. S. Dean and A. 
Lefevre, Phys. Rev. Lett. 86, 5639 (2001). H. A. Makse 
and J. Kurchan, Nature 415, 614 (2002). A. Fierro, 
M. Nicodemi and A. Coniglio, Europhys. Lett. 59, 642 
(2002). G. De Smedt, C. Godreche, J.M. Luck, Eur. Phys. 
J. B 32, 215-225 (2003). 

[15] M. Nicodemi, A. Fierro, A. Coniglio, Europhys. Lett. 
60, 684 (2002). A. Fierro, M. Nicodemi and A. Coniglio, 
Phys. Rev. E 66, 061301 (2002). A. Fierro, M. Nicodemi, 
A. Coniglio, Jour. Phys.: Cond. Matt. 15, S1095 (2003). 

[16] M. Dijkstra, R. van Roij and R. Evans, Phys. Rev. E 59, 
5744 (1999). M. Fasolo and P. Sollich, Phys. Rev. Lett. 
91, 068301 (2003). 

[17] M. Bose, P.R. Nott and V. Kumaran, Europhys. Lett. 
68, 508 (2004). An effective attraction between two large 
grains was also found in a 2d vertically vibrated granular 
bed of smaller ones in Ref. DA. Sanders, M.R. Swift, 
R.M. Bowley, and P.J. King, Phys. Rev. Lett. 93, 208002 
(2004). 

[18] For a given grain configuration, {nf}, n({nf}) 
has a tractable expression: n({nf}) = 

lim/f^oo exp {-KHEdw} where HEdw 



y. 



i n 6 z 

n . 



[19] M. Tarzia, A. de Candia, A. Fierro, M. Nicodemi and A. 

Coniglio, Europhys. Lett. 66, 531 (2004). 
[20] A. Fierro, M. Nicodemi, M. Tarzia, A. de Candia and A. 

Coniglio, cond-mat/ '0412120. 
[21] M. Mezard and G. Parisi, Eur. Phys. J. B 20, 217 (2001). 

G. Biroli and M. Mezard, Phys. Rev. Lett. 88, 025501 

(2002). 

[22] Along with the fluid solution, also a crystalline solution 
characterized by the breakdown of the horizontal trans- 
lational invariance, can be found. Here we consider the 
case where the crystallization is avoided. 

[23] /C&° = A^/S^ where A ( n z) = 1 + 8 n ,ie s i ) and S (z) 



Moreover we have that P„ 



A[ z 

5„, 2 e"i ' +e v 2 ■ , R\ z > = 5 n , ie a i ' + <5„, 2 e c i ' + e c 2 ' , Q 



1 + S n . 1 e u i ) + 

(*) „(*) 



l + e a i [l+e u 2 "} + e d 2 " [l+ e T "' +e u 2 ' "]+R% 
[24] The free energy, F, is given by F = J2f=o Ai? " (z) - 
(fc-i) A r(z) v-ff-iAzr« _.v. „-af« 

(5 ,(z) ) fc ~ 1 Q (z) 



M 

',2 
+ 



T,z=o AF i,i> where - e 

7< Z > „ (z) + + 

= 1 + E„=i, 2 e c " (l + * n , 2 e"i +e4 K 

(z+l) 



.-<*) 



= l + 2e s i +2e s 2 + e z 



) + 
(«) 



[25] Each tap consists in raising the bath temperature to a 
finite value and, after a lapse of time quenching the bath 
temperature back to zero. The measurements are per- 
formed when the shake is off and the system is at rest. 



